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For the model of frustrated spin-1/2 Heisenberg magnet described in A. A. Nersesyan and A. 
M. Tsvelik, (Phys. Rev. B67, 024422 (2003)) we calculate correlation functions of staggered mag- 
netization and dimerization. The model is formulated as a collection of antiferromagnetic chains 
weakly coupled by a frustrated exchange interaction. The calculation done for the case of four 
chains demonstrates that these functions do not vanish. Since the correlation functions in question 
factorize into a product of correlation functions of spinon creation and annihilation operators, this 
constitutes a proof that spinons in this model propagate in the direction perpendicular to the chains. 

I. INTRODUCTION 

In Q Alexander Nersesyan and one of the authors presented a proof of existence of fractional quantum number 
excitations in a model describing a certain frustrated magnet in the number of dimensions greater than one. This 
magnet consists of spin-1/2 antiferromagnetic Heisenberg chains weakly coupled by a frustrated antiferromagnetic 
interaction: 

j,7l ' Sj_J_X,7l 

j,n I fi=±l ) 

where Sj in are spin-1/2 operators, and Jn >> J r , Jd > 0. The interaction pattern of this model resembles the flag 
of American Confederation which gave the model its name. Fractional quantum number excitations appear when 
J r (n,n + fi) = 2J ( i{n,n + /i). If this condition is fulfilled, the interaction between staggered components of the 
magnetization on chains n and n + /i vanishes. The weakness of the interchain coupling allows us to employ the 
continuous description. In the continuum limit each Heisenberg chain is represented by the SUi(2) Wess-Zumino- 
Novikov-Witten (WZNW) model and the relevant part of the interchain interaction is reduced to the interaction of 
spin currents with different chirality: 



JV 



n=l 



H n + 7(71, n + fi) J da;(J n J n+(U + J n+M J r . 



(2) 



Here H n is the SUi(2) WZNW Hamiltonian 



H n = ^Jdx(:3l: + :3l--) (3) 

where v = irJ\\ao/2 is the spin velocity, J%(x), J%(x) are operators representing holomorphic and antiholomorphic 
currents belonging to the SU(2)i Kac-Moody algebra. The currents with different n commute. The coupling constant 
7(n, n + (i) ~ J r (n, n + //). Continuum limit Hamiltonian J2J coincides with the one introduced by Emery, Kivelson 
and ZacharQ in the context of theory of stripes. We emphasise that the interchain interaction in this frustrated model 
remains relevant and generates spectral gaps. Therefore this model does not belong to the class of sliding Luttinger 
liquid models where soliton excitations remain confined to the chains. 

Though the authors of [l| have managed to proof the existence of spin S — 1/2 excitations in model ©, it remained 
unclear whether these excitations are able to propagate in the direction transverse to the chains or remain confined. 
In the Confederate Flag model the staggered magnetization operator creates two non-interacting spinons (the absence 
of interaction follows from the fact that these spinons belong to different sectors of the Hamiltonian H + , H ~ , -see 
below) . Therefore one way to resolve the problem of propagation is to calculate correlation functions of the staggered 
magnetizations between different chains. If these correlation functions do not vanish, the excitations propagate. The 
difficulty is that such correlation functions are essentially non-perturbative objects (they remain zero in any order of 
perturbation theory) and to calculate them one has to somehow go beyond perturbation theory. The corresponding 
methods are available for a finite number of chains (namely N = 2,3,4 where exact solutions are available). Though 
the most interesting case corresponds to an infinite number of chains N — > oo, information obtained for finite may 
also provide valuable insights. The calculations of the interchain correlation function of staggered magnetizations done 
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in for the case of two chains gave a nonvanishing answer, but one may argue that N — 2 is too small a number 
to allow even a qualitative extrapolation to the thermodynamic limit. Similar calculations done for N — 4 could not 
distinguish between intra and interchain correlation functions. In this paper we perform accurate calculations for the 
case of four chains N = 4. We find that the interchain correlation functions do not vanish. 

In the continuum limit the staggered magnetization N = (— l) J Sj and the dimcrization e = (— 1) J (SjSj+i) operators 
become smooth fields. For the chain number n they are expressed through matrix elements of the 5 = 1/2 primary 
field of the n-th WZNW model: 

.g CTC r'(n; x) = <We(«; x) + i(o% CT /N(rc; a;) (4) 

The action for model J2J is local in g and posseses the global SU.r(2)xSUl(2) symmetry. Namely, the action remains 
invariant under the following transformations: 

92n — ► 92nV, ,92m+l — * V + g 2m +l> 

9ln — * Ug2n, 32m+l ~> 92m+lU + , (5) 

where V, U are coordinate-independent SU(2) matrices. This symmetry dictates the following form of the two-point 
correlation functions: 

((9 CT ia 2 (r,£;2n)p+ CT4 (0,0;2?7i))) = S sia J 82tT3 V 2nt2m (T, x) 

({9n<T 2 (r,x;2n)g a3 o- 4 (0j 0; 2m + 1))) — 8 sirT4 6 S2(7:j 'D2n,2m+l ( T > x) (6) 

where T> do not contain spin indices. Substituting Eq.(0J) in the above equations we obtain the following important 
relation between correlation functions of staggered energy density and staggered magnetization: 

({e2n(r,x)e 2m (0,Q))) = (Wn(^)C(M))) 
((e 2 „+i(r I a:)e 2m+1 (0,0))) = ((N^ +1 (r, x)N% m+1 (0, 0))) 
((e 2n (r,x)e 2m+1 (0,0))) = -((N« n (r, x)N% m+1 (0, 0))) 

(7) 

Now let us recall the important feature of Hamiltonian J2J): as was noticed in 0], it separates into a sum of two 
commuting parts: H = H + + H~ (such structure of the Hamiltonian was first observed in the context of Kondo 
lattice in |3) ■ Therefore eigenstates are separated into two sectors with different parity. The Hamiltonian density in 
the plus parity sector is 

n + = ^ (: Ji 2 : + : J 3 2 : + : J 2 2 : + : 3\ :) + (AJi + A'J 3 )(AJ 2 + A'J 4 ) (8) 

For the purposes of this paper we find it convenient to consider different copling constants A / A' such that the 
difference between these coupling constants is small (|A — A'| << A). Physically this means that we still consider the 
frustrated interchain interactions, but allow their amplitudes to vary between different chains. 

By using the Abelian bosonization procedure for the SUi(2) currents one can write the Hamiltonians H in terms 
of bosonic fields (p, (p living on each chain (see for the details). In this representation fields tp\$ (<?i 3) interact with 
the fields (p2A (^2,4) such that 

H = H+[Lp 1}3 ; if 2,4} +-ff~[£i,3;V2,4] ( 9 ) 
The convenience of this representation is related to the fact that it allows to represent g(j) in a factorized form: 

1 



~C aa t : expf-iv^Mcr'/ 3 + cr'ip)] := C, 
v2 



(7(7' ^<T^(T' j 



C aa , = eM 1 "™')/^ (10) 

where 

Zcr = exp[iCTv27r</'], z CT = exp[— 'kt^/2tt0\, (a = ±1). (11) 
Eas. H9llQ(l select a set of potentially non- vanishing interchain correlators: 

G13 = (( e i ^(i) e -i^(3))) ) Q 13 = (^ e i^(i) e -i^(3))) 

^12 = ((e i ^(i) e i^(2))) ; d 12 = (( e i^(i) e -i^(2)^ (12 ) 
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From these correlators one can obtain correlation functions of both vector N and scalar e staggered fields: 
(<e 2n (T,z)e 2m (0,0))) = ((N$ n (T,x)N2 m (0,0))) = |G 2 „, 2m (r, x)\ 2 

((e 2 „(r,x)e 2m+1 (0,0)» = -«JV&(r, x)N$ m+1 (0, 0))) = \D 2n , 2m+1 (r, x)\ 2 (13) 

According to these formulas spins on neighboring chains are oriented antifcrromagnetically. We emphasise that 
this conclusion follows just from the symmetry arguments and the chiral decoupling of the operators (|1C)|> and the 
Hamiltonian ijjj}. 

II. EXACT SOLUTION OF MODEL © 

The exact solution of model JSJ) was obtained in The spectrum includes heavy solitons and antisolitons with 
mass M and a light singlet Majorana fermion with mass m. The mass ratio is m/M ~ (A — A') and vanishes in the limit 
of equal couplings. This is the limit which was studied in 0, it corresponds to the uniform interchain interactions and 
periodic boundary conditions in the transverse direction. The presence of such periodicity constitutes an additional 
symmetry whose presence leads to vanishing of certain correlation functions. Therefore we prefer to keep m/M finite. 
The two-particle scattering matrix of model (JSJ) is given by 

V ° a >" e"/ 2 +i 1 / 

where 5[su(2)] is a 2x2 scattering matrix of the SU(2) Thirring model. 
The Thermodynamic Bethe Ansatz (TBA) are 

Tm f TM f 
F/L = / d6»cosh<91n[l + c ei W/ T ] / d0cosh01n[l + C ^/ T ] (15) 

2?r J 27T J 

e n (9) = T S *ln[l + e e "- 1 W/ T ][l + e £ "+ 1 W/ T ] 

— 5 n ,i m cosh 9 — 8 n ^M cosh 9 (16) 

These equations are valid for A' w A. 

It is worth mentioning that model JSJ is generated in the relativistic limit as the spin sector of the following fermionic 
model: 



H = dx 



2 

a = ±1/2,3 = 1,2; 5 >0 (17) 
From the solution of this model obtained in |5| it can be extracted that 



m 7r / tj_ 
M ~ 8 V^F 



2 



\n 2 (e F /M) 



TV 

which provides a direct relationship with model JSJ 



M = ^kpgc^i-rrkp/g) (18) 



III. THE CORRELATION FUNCTIONS 



We shall calculate the correlation functions using the formfactor approach. In this approach one calculates matrix 
elements (formfactors) of various operators using their transformation properties under various symmetries of the 
problem. General information about the method can be obtained from various review articles @,|3- 
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Let us consider operators exp[ier\/27n^i] and exp[i<7v27ry>3]- We will be interested only in the limit m << M, 
therefore it is sufficient to consider the states with one soliton (antisoliton) and arbitrary number of light particles. 
In this case it is reasonable to suggest the following formfactor expansions: 



cxp[ierV27npi]|0 >= Q even + Qodd 

ex.p[iaV2n<P3]\0 >= Q even - Q odd (19) 

where Qeven (Qodd) has matrix elements with the states with of one heavy soliton and even (odd) number of light 
Majorana fermions. Notice that the on-chain correlation functions for both operators are the same as it must be. 
The problem is simplified by the fact that the Majorana fermions have a non-trivial scattering matrix only with the 
heavy particle. Otherwise the dependence of the matrix element on the Majorana fermion rapidities is like for the 
Ising model. Using the results of (g, where a systems with a similar S-matrix was studied, we obtain for the following 
expressions: 

2n 

< s(fa; X i(Pi), -X2n(02n)\Q even = (C+^V/ 4 J] ^((3 - fa) JJtanh (/%/2) (20) 

J=l i-3 
2n+l 

< S (f3); X l(fa),-X2n + l(fan + l)\Qodd = CT V V/ 4 J] V^H/3-/^II tanh (/V2) (21) 

J=l i-3 

where (3 is the soliton's rapidity and fa stand for the rapidities of the light particles. C ± are normalization factors to 
be determined later. The function ip(~\fa satisfies the following equation 

^-)^i-){P-m) = -l_ (22) 



and is given by 



^~\I3) =e- /3/4 V (0) (/3), (23) 

^)(^) = 2-^ex P {- r^2sin^p + i,W2] + sinh 2 W2) 
y J uj 2 sinh(7rw) cosh(7rw/2) 

As follows from (|22|l . this function has the following asymptotics: 

^\/3^+oo) ie"^ 2 

V> (_) (/?^ -oo) -> 1 (24) 

The operators exp[iv27T¥>(2,4)] are expressed in a similar way, but with all rapidities having the opposite sign. 
Therefore their expansions include the function ijj^ + \f3) = ip^'{—fa) which decays at — oo and approaches 1 at 

Let us briefly comment on expressions H2UI21|) . The operators under consideration have Lorentz spin 1/4. This 
dictates that under the Lorentz transformation /3 — > /3 + 9, f3j — > fa + 6 the matrix element of that operator must 
acquire a factor exp(0/4). As we see from (|20I21|I . this property is fulfilled. Another property is that the formfactor is 
multiplied on the S-matrix S((3 7 fa) whenever two particles with rapidities f3 and fa are interchanged. Since tanh(/3y/2) 
is an odd function, this is obviously fulfilled for the Majorana fermions. The function -0 plays the same role for the 
ineterchange of a Majorana fermion and the heavy particle. The choice of asymptotics of functions tp^ is determined 
by the fact that in on small distances bosonic exponents (|19J) become (anti)holomorphic fields. 

Substituting l|2UI21|l into Eas. l|19|) and performing certain algebraic manipulations (see Appendix 1), we arrive at 
the following expressions for the correlation functions on the same and on different chains (r 2 = r 2 + x 2 , we put 
v = l): 

G pq = ((eiv^Mrie-iv^M)) = ( 1— ^ C(M i m) x/% J : pq (r) 

\T + IX J 

D pq ee (( e i^W,jiV5^(!))} = C(M 4 m) 1/8 ^ pg (r) (25) 
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where C is a numerical coefficient, p, q = 1, 3 for G and p = 1, 3; q = 2, 4 for D. 
_Fu(r) = ^33 (r) = 

= / d0eP/ 3 e- Mrcaih(, exp | £ - / d/3 fj |V ( - } (/3 - /3 J )| 2 e- mrcosh * J|[cosh(A J /2)]" 1 1 (26) 
^13 (r) = ^24 (r) = 

= / d/3e^ 2 e- Mrcosh/3 exp J ^ / d/8 JJ (/? - /3 J )| 2 e- mrcosh * JJ[cosh(/3^ /2)]" 1 } (27! 

For the correlation functions on neighboring chains (sec Eas. (|13|l (p — 1, 3; q — 2, 4) we have 



D 12 (r) = / d/3e- MrcoBh ^exp| £ ± f d/?f[ |^ (0) (/3 - &)| 2 e- mrcoshft ^coshC^)]" 1 
Du(r) = y d/3e~ Mrcosh ^exp j ^ / ^ft I^G 3 " M\ 2 e~ mrcoBh ^ Y\[cosh(f3 tJ /2)}- 1 



(28) 



n=l ' j i>j 



The normalization factor (M 3 m) 1//S in Eas. l|25|) is dictated by the facts that (i) the single chain correlation functions 
do not vanish in the limit m — > 0, (ii) the scaling dimension of the operator is 1/4. In the limit Mr >> 1, mr << 1 
the integral in (3 converges at |/3| ~ (Mr) -1 / 2 . Since |V* (AOI 2 ~ ex P( — I/?? 1/2), the integrals in /3j converge even at 
mr = 0. As result we get 

£>ia(r) = C(Af 3 m) 1 / 8 X (Mr), D 14 (r) = C(A/ 3 to) 1 / 8 ^ (M^) (29) 

where is C the numerical constant which remains undertermined. To calculate the other asymptotics is a more 
complicated task. The calculations are done in Appendix 2; the result is 

^• u (r)=C u M- 1 /V Mr r- 5 /8 

F 13 (r) = Cism^M-Wr-We-*"- (30) 

where C are again undertermined numerical constants. From here we derive the following estimates for the singularities 
in the imaginary part of the staggered magnetic susceptibility (dimerization) at s = 2M, (s 2 = lu 2 — q 2 ): 

9fmxu(«) ~ (s- 2M)- 1/4 , 3mxi 2 m 1/4 {s ~ 2M)~ 1/2 , 3mxi 3 (s) - mAf" 5/4 6»(s - 2M) (31) 

These expressions are valid at s — 2M >> m so that 3mxi2 cannot really become greater than Smxu. 



IV. CONCLUSIONS 



Looking at expressions for the interchain correlation functions, we see that thay are proportional to powers of the 
particle masses. These masses are generated dynamically in the theory, they are exponentially small in the inverse 
coupling constant and therefore the interchain tunneling of the solitons is an essentially non-perturbative process. This 
conclusion coincides with the result for two chains ,lj. The four chain case, however, introduces a new feature: the 
presence of the singlet particle. It is interesting that the prefactor in the interchain correlation functions contains both 
masses and vanishes when the mass of the singlet particle goes to zero (this occurs for periodic boundary conditions 
in the transverse direction; the case considered in This indicates that the singlet sector plays an important 

role in the interchain propagation of solitons probably providing conditions for an uninhibited tunneling. We remind 
the reader that the fraction of the Hilbert space occupied by singlet excitations grows with the number of chains pj 
and these excitations will certainly play an important role in the thermodynamic limit determining interchain soliton 
propagation. 

In conclusion we would like to emphasise a curious symmetry between correlations functions of the staggered energy 
density and magnetization (|13f) existing in Confederate Flag magnet. This symmetry unites e and N into a four- 
component vector giving rise to an analogy between this model and models of non-collinear magnets. As it was 
pointed out in [3, spinons naturally exist in disordered non-collinear magnets and it looks likely that Confederate 
Flag model provides a microscopic realization of that scenario. 
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APPENDIX 1 



G n = ^2— } / d/3e /3/2 e- MTCOsh/3 - iM;i:sinh/3 Y[d(3j [J - /^l V mTCOBh *- ixmsillh * J^tanh 2 (/3 y -/2) 

ra=0 ' j j=l i<j 

Gi 3 = J2^~ / d/3e^ /2 e- MTCOsh ^- 1Ma:sinh ^JJd^ JJl^-^/?- / 3 J )| 2 e- mTCOsh/3 ^ lxmsinh ^ []tanh 2 (/3y/2) 



n=0 j j=l i<j 



By the shift of the integration contours one obtains the following expressions: 

1/4 



(32) 



G pq (r, x) = — — F„(r), t' = t z + X* (33) 

\T + IX J 

^ii(r) =^ 33 (r-) = 

^ P n n 

12 —} / d/3n d ^ e/3/2e " Mrcosh/3 n - ^-)i 2 e" mrcoshft n tanh2 (^/ 2 ) 

n=0 ' j j=l i<j 

= f d^ 2 e~ Mrcosh ^exp \ - I df3f[\^-\p - (3 J )\ 2 e- mrcosh ^Y[[cosh(l3 l] /2)}- 1 I (34) 
J (n=i nJ j i>j J 

^13 (r) =^ 3 i(r) = 

/ d^dfte^e- 1 '" 08 ^ J] - &)| 2 e- mrcosh ^ [] tanh 2 (/%/2) 



n=0 " j j = l i<j 

= j d@eV 2 e- Mrcosh Pexplj2^- J dpf[\^- ) (p-(3 3 )\ 2 e~ mrcosh ^Y[[cosh((3 lJ /2)}- 1 



(35) 



Here we have used the fact that the correlation functions are equal to the determinants of the Fredholm operators 
(see Appendix 2) and used Ea, 137|l . The correlation functions on neighboring chains are given by similar integrals: 

D pq = ((eiv^Mei^^ 2 ))) = ^( r ) ( 36 ) 

i /• n n 

^i 2 (r) = ^2 nil d / 3 Il d ^ e " MrCOSh/3 II ~ ~ /3,)e- mrcosh * n tanh2 (^/2) 

n=0 ' j j=l i<j 

I -i \n f n n 

= E —T j d/3n d ^ e " MrC ° Sh/3 II - - /3,)e— h ^ \{t^ 2 {f3 l3 /2) 

n=0 j j=l i<j 

Taking into account that ip + (P)^^ (/?) = |V> (0) (P)\ 2 and using Eq.^ we arrive at Eas. (|28l) . 
APPENDIX 2 

In this Appendix we calculate the asymptotics of correlation functions T pq {r). First, let us consider as an example 
the scaling of the Ising model. It is well-known that the correlation functions of the order and disorder parameter 
operators for this model can be expressed as follows: 

(fj,(x)(j.(0)) = G+(mr) + G_(mr), 
(a(x)a(O)) = G+(mr) - GL(rw), 
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G±(p) = ^(±r^ / dA--- / d^ n e-^ cosh ft J] tanh 2 ^(A -ft) 

"=° -oo -oo «i 



The functions G± are determinants of the Fredholm operators: 

G±(p) =det (/ ± K) 
where K is the integral operator with the kernel: 



cosh - /3 2 ) 



From here one can conclude that G± satisfy the Painleve equation ^} and perform a rather detailed analysis of these 
functions. The analysis of asymptotic behavior is a simplier task however and can be performed as follows. Using the 
formula 



log(detA) = Tr(logA) 



(37) 



one obtains: 



oo i r p n 

log(G±(p)) = ^(±r- / dA-"- / d^e-^coshft-Q 

« =1 oo 

where /3 n +i = (3\. Consider one of integrals: 

oo oo n 

f djSi--- / d/3„e-^ cosh ft -Q . 



^ cosh 2 |(/3j - A+i) 



cosh 2 i(/3j - 

We are interested in the asymptotics for p — > where the integral diverges. Let us introduce new variables: 



(38) 



Ai = Pi-0 i+ i, i=l,-",n-l, 

N 

ui = ^ cosh(/3j) 

3=1 



(39) 



Then the integral becomes 

oo 

/dA, 

— oo 

where 



dA n _i 2 / dw 



n 



D(A) 



y/u* - D(Ay " cosh 



cosh 



D(A 1 ,... > A„_ 1 )=(^eft')(^e- 



It is clear that the integrals over A.,- are always rapidly converging while the integral over u> diverges at p — ► 0. The 
estimation of leading contribution is straightforward. We give the final result for G + which diverges at p — > 0: 



MG+(P)) 



oo » n. n— 1 

2 ^/ da ,.../ da „_ in 

3 = 1 -oo -oo 3 = 1 COSn I 2 



1 



i cosh (2) cosn ( ^2^ ) 



log(p) 



(40) 



The expression in square brackets is the anomalous dimension. It equals \. This value is known from many sources, 
but actually even direct summation of series is possible. In a similar way we obtain 



log(G_(p))~p 3 / 4 



(41) 
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Now let us turn to the case considered in the present paper. 



oo n 

E(±) n i /^i- / d J a B e-^=^ftJJtanh a i(A- J 8 i )n|^ ( - ) (^i) 
"=° '-"L -"oo *<i i =1 



The functions .F± are determinants of the integral operators 

T±(p) = det (J±L) 

where L is the integral operator with the kernel: 



_ „— jp(cosh/3 1 +cosh / 3 2 )_ 



1 



cosh - (3 2 ) 



So, similarly to the Ising case the logarithms of these determinants can be written as 
log(^±(p)) = 

„ OO OO 

1 . . 



p cosh j3 : 



n 



\ cosh 2 | (ft - 



Consider the integral 



d/3 



l • • • 



d/3„e' 



- p ^ cosh /3j 



n 



\ cosh 2 |($ - ft+i) 



and make change of variables l|39[) . The difference with the Ising case is due to the fact that we need to express f3j in 
terms of the new variables and to substitute them into tp. We have: 



P 3 = log (u T y/u; 2 -D 2 {A)^j - log (e~^ ^ e 0k ^j 



(42) 



where the last term depends only on A's. The integral over u> is taken between the limits -D(A) and oo over two 
branches: one the first one one takes + sign in l|42(l . and on the second one the minus sign. The part divergent at 
p — * comes from the region where — log(w) is large. It is not the second branch of the integral over u>, since |?M - )(/3)| 
is rapidly decreasing as (3 — > +oo. So, only the first branch contributes, and, since 

I^G 9 )! rapidly approaches 1 as 

(3 — > — oo, we can replace all by 1. The result of this considerations is that the asymptotics of the integral 

is \ of what we had in the Ising case. So, 



log (F + (p)) ~ - i log(p), log (JF_( p )) ~ | log(p) 



(43) 
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